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We introduce a new expression for the Z2 topological invariant of band insulators using non- 
Abelian Berry's connection. Our expression can identify the topological nature of a general band 
insulator without any of the gauge fixing problems that plague the concrete implementation of 
previous invariants. The new expression can be derived from the "partner switching" of the Wannier 
function center during time reversal pumping and is thus equivalent to the Z2 topological invariant 
proposed by Kane and Mele. 
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I. INTRODUCTION 

Topological invariants play a very important role in the 
classification of band insulators. The studies on integer 
quantum Hall effect (IQHE) show that 2D band insula- 
tors without time reversal symmetry can be classified by 
the Chern number - an integer describing the topological 
structure of a set of fully occupied Bloch bands with- 
out Kramers degeneracy. Systems with non-zero Chern 
number exhibit IQHBii^. 

A similar idea can be also applied to band insulators 
with time reversal symmetry. Recently a Z2 topologi- 
cal invariant has been proposed by Kane and Mele to 
characterize the time reversal invariant band insulators 
in 21>2.ii. According to this new topological invariant, all 
the 2D band insulators with time-reversal invariance can 
be divided into two classes. The normal insulators with 
even Z2 number and topological insulators with odd Z2 
number—"—. The 2D topological insulators will exhibit 
a quantum spin Hall effect (QSHE)^i^, which is charac- 
terized by the presence of helical edge states^"—. In- 
terestingly the Z2 topological invariant can also be gen- 
eralized to the 3D band insulators with time reversal 
symmetry^i^ii^. In this case, there are four independent 
Z2 topological numbers: one strong topological index and 
three weak topological indicesi^"— . The 3D time rever- 
sal invariant band insulators can be classified as normal 
insulators, weak topological insulators (WTI) and strong 
topological insulators (STI) according to the values of 
these four Z2 topological indices. Among them, the STI 
attracts much attention due to its unique Dirac type sur- 
face states and robustness against disorder—"—. The he- 
lical spin structure of the Dirac type surface states has 
been experimentally implied by the standing wave struc- 
ture in STM images around an impurity scattering cen- 
ter and measured directly by spin resolved angle resolved 
photo emission spectra (ARPES)^"— . The Dirac type 
2D electron gas living on the surface of STI or at the 
interface between STI and normal insulators provides a 
new playground for spintronics and quantum computing. 



Since the Z2 invariant characterizes whether a system 
is topologically trivial or nontrivial, its computation is 
essential to the field of topological insulators. For band 
insulators with extra spacial inversion symmetry, the Z2 
topological numbers can be easily computed as the prod- 
uct of half of the parity (Kramers pairs have identical 
parities) numbers for all the occupied states at the high 
symmetry points^^. The situation becomes complicated 
in the general case where spacial inversion symmetry is 
absent. At the present, numerically there are three dif- 
ferent ways to judge if a band insulator without inver- 
sion symmetry is a TI or not. i) Compute the Z2 num- 
bers using the integration of both Berry's connection and 
curvature over half of the Brillouin Zone (BZ). In order 
to do so, one has to set up a mesh in the k-space and 
calculate the corresponding quantities on the lattice ver- 
sion of the proble m^^i^^ . Since the calculation involves 
the Berry's connection, one has to numerically fix the 
gauge on the half BZ, which is not easy for the realis- 
tic wave functions obtained by first principle calculation, 
ii) Start from an artificial system with spacial inversion 
symmetry, and then smoothly "deform" the Hamiltonian 
towards the realistic one without inversion symmetry. If 
the energy gap never closes at any points in the BZ dur- 
ing the "deformation" process, the realistic system must 
share the same topological nature with the initial ref- 
erence system whose Z2 number can be easily counted 
by the inversion eigenvalue formula. Unfortunately mak- 
ing sure that the energy gap remains open on the whole 
BZ is very difficult numerically, especially in 3D. iii) Di- 
rectly calculate the surface states. For most of the TI 
materials, the first principle calculation for the surface 
states is numerically heavy. Therefore it is valuable to 
develop a mathematically equivalent way to calculate the 
Z2 numbers of a band insulator, which satisfies the fol- 
lowing conditions: first it should use only the periodic 
bulk system; second, it should not require any gauge fix- 
ing condition - thereby greatly simplifying the calcula- 
tion; third, it should be easily applied to general systems 
lacking spacial inversion symmetry. 
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In the present paper we propose a new equivalent ex- 
pression for the Z2 topological invariant using the U(2N) 
non-Abelian Berry connection. Based on this new ex- 
pression, we further propose a new numerical method to 
calculate the Z2 topological number for general band in- 
sulators, without choosing a gauge fixing condition. The 
main idea of the method is to calculate the evolution 
of the Wannier function center directly during a "time 
reversal pumping" process, which is a analog to the 
charge polarization^!^. We derive that the center of 
the Wannier function for the effective ID system can 
be expressed as the U(l) phase of eigenvalues of a ma- 
trix obtained as the product of the U(2N) Berry connec- 
tion along the "Wilson loop". The Z2 topological num- 
bers can be expressed as the number of times mod 2 of 
the partner switching of these phases during a complete 
period of the "time reversal pumping" process. Using 
this new method, we have recalculated the Z2 topolog- 
ical numbers for several TI systems, including strained 
HgTe, Bi, Sb and Bi2Se3, and found the "partner switch- 
ing" patterns, which differentiate between topologically 
trivial and nontrivial behavior. The rest of paper will 
be organized as follows: in section II we derive the new 
mathematical form of the Z2 numbers through the "Wil- 
son loop"; we apply the new method to various TI sys- 
tems in section III; we prove the equivalence of the new 
methods and the Z2 number propose by Fu and Kanei 
in the appendices. 

II. THE FORMALISM 

The Bloch wave functions describing the band struc- 
ture of a translationally invariant system can be ex- 
pressed in terms of a complete set of local basis labeled 
by the unit cell i and some other quantum number a as: 

|^„fc(r)) = ^Y.u^^{VY^-^^ |</'a(r-R,)) (1) 



N 

We will first focus on the 2D system. The topological 
nature of a 3D insulator can be determined by looking at 
these effective 2D systems with fixed hi = and hi = tt 
(i = x,y, z). 

Following Fu and Kane, for a 2D insulator we can 
further reduce the dimension by fixing ky and study 
the "time reversal" pumping process with the adiabatic 
change of ky. The main idea of our new formalism is to 
directly look at the evolution of Wannier centers for these 
effective ID systems in the subspace which contains only 
occupied states. Fixing ky, we can define the position 
operator for the effective ID system as 

X = J2 e'''^-^'|</>a(r - R.))(0c(r - R.)| (2) 

where 6kx = j^^, is the number of real-space unit 
cells along the x Jirection, ax is the lattice constant, a is 



the orbital and spin index and R; labels the unit cell. The 
operator A is a unitary operator with all the eigenvalues 
being g''^'^'"'^', whose phase represents the position. The 
eigenvalue of the position operator can be viewed as the 
center of maximum localized Wannier function (MLWF) 
formed by the bands included in the operator X. Be- 
cause the local basis set a is assumed to be complete , 
such MLWFs are always well defined. As pointed out by 
Fu and Kane, the Z2 topological invariant can be deter- 
mined by looking at the evolution of the Wannier func- 
tion center for the effective ID system with fixed ky in the 
subspace spanned by the occupied bands only . There- 
fore we should consider the eigenvalue of the projected 
position operator defined as 



Xp = PXP 



nmgo k^k'^,ia 
X |0a(r - R,))((/)a(r - R.O|V'mfeifcJ(V'mfcifeJ 

nm^o kj-k'A 



E/ ^na{^x)Uma{kx) 

5{kx + 5kx - k'^) ^ \i'nk^ky){lpr, 



(3) 



where o means the occupied bands. The above operator 
can be written in a more suggestive matrix form 



Xp{ky) = 
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(4) 

where F^^^{ky) = Y,a'^na{kx,i,ky)uraa{kx,i+i,ky) are 
the 2M X 2M matrices spanned in 2M occupied states 
and kx^i — are the discrete k points taken along the 
X-axis. 

The eigen problem of Xp{ky) can be solved by the 
transfer matrix method. We can define a product of 
Fi^i+i as 



D{ky) — i^0,l-P'l,2^2,3 • • • ^Ar^-2,7Vx-l-^iVx-l,0 (5) 

D{ky) is a 2M x 2M matrix, which has 2M eigenvalues: 



m= 1,2,--- ,2M 
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where 9^ is the phase of the eigenvalues: 



A" 



(6) 



We can easily prove that the eigenvalue of projected 
position operator Xp can be simply related to the eigen- 
value of the above D-matrix by 



Af 



(7) 



where n ^ 1, 2, • • • , Nx- We can further prove that the 
D-matrix is unitary and all the |A^| equals one. 

The evolution of the Wannier function center for the 
effective ID system with ky can be easily obtained by 
looking at the phase factor 0^. At fcj, = 0, the eigen- 
values of the D-matrix appear in degenerate pairs due to 
time reversal symmetry, which results in pairs of Wan- 
nier centers sitting at fcj, = 0. When ky moves away from 
the origin, the Wannier center pairs split and recombine 
at ky = TT, as shown in FigdJ Because 6^ is a phase 
factor, when two 6'^s meet together, they may differ by 
integer times of 27r. Therefore the evolution of each Wan- 
nier center pair will enclose the whole cylinder an integer 
times, which can be viewed as the winding number of 
the Wannier center pair. The Z2 topological number is 
related to the summation of the winding numbers for all 
the pairs. If it is odd, then the Z2 topological number 
is odd. It seems that the total winding number of the 
Wannier center pairs should generate an integer class Z 
instead of Z2. To clarify this point, let's look at the evo- 
lution of a Wannier center pair with winding number An. 
In that particular shown in FiglljC), the pair 

of Wannier centers must have an extra "accidental" de- 
generacy between ky = and tt, which is not protected 
by any symmetry and can be removed by "deforming" 
the Hamiltonian slightly to make the crossing of the lev- 
els become an anti-crossing as shown in FiglljC). The 
"deformation" process will thus change the total wind- 
ing number by Air and make it 0. Therefore only the total 
winding number mod 2 is a topological invariant. 

The eq.([5|) can be viewed as the discrete expression of 
the Wilson loop for the U(2M) non-Abelian Berry's con- 
nection. It is obviously invariant under the U(2M) gauge 
transformation and thus can be calculated directly from 
the wave functions obtained by first principle method 
without choosing any gauge fixing condition, which is 
the biggest advantage of the present form of the Z2 in- 
variance. The equivalence to the Z2 number proposed by 
Kane and Mele will be proved rigorously in the appendix. 



III. NUMERICAL RESULTS 

In the present section, we will implement the method 
and explicitly compute the Z2 invariant for a series of 
systems. For each particular system, we calculate the 
evolution of the 9^ defined in equation® as the function 



Ay 




Figure 1: Schematic plots of the Wannier center curves: (A) 
for the trivial case, the Wannier center winding the cylinder 
zero times; (B) the Wannier center winding the cylinder one 
times; (C) the Wannier center winding the cylinder twice and 
the cross point labeled by P is not protected by time reversal 
symmetry and it is usually eliminated by some perturbation 
terms. (D)the Wannier center winding the cylinder 3 times, 
which is topologically equal to the case in B. 



of ky from zero to tt. The winding number of the Wannier 
center pairs defined in the above section can be checked 
in an equivalent way which is much simple in practice. 
We first draw an arbitrary reference line parallel to the 
ky axis, then compute the Z2 number by counting how 
many times the evolution lines of the Wannier centers 
crosses the reference line. 



A. BHZ model 

Bernevig Hughes and Zhang (BHZ) showed that for 
an appropriate range of well thickness, the HgTe/CdTe 
quantum well exhibits an inverted sub-band structure. In 
this inverted regime, the system exhibits a 2D quantum 
spin Hall effect^i^. BHZ introduce a simple four band 
tight binding model to describe this effect: 



Heff{kx, ky) 



H{k) 






i/*(-k) 



(8) 
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where i?(k) 
i sin k,,],d 



= £(k) 

-2B[2- 
- cos k,. 



- ^[sinfc^ 
.,j,e(k) = C 



(ii(k)(Tj, di + idi 

j,j ,1^3 = — yb ~ '^'^^ ~ '^'-'^ 
2Z)[2— cos /c^c — cos ky]. Real HgTe does not have inversion 
symmetry but the BHZ toy model does. To describe the 
inversion symmetry breaking effect we add a new term 



H' 



A 

-A 

-A 

A 



(9) 



We apply the new method to calculate the shift of 
the Wannier function center based on the above model 
Hamiltonian and show the results in Figl2] 




Figure 2: Wannier centers for the BHZ model. (A) For 
the QSH phase (A=-13.68, B=-16.9, C=-0.0263, D=-0.514, 
M=-2.058, A=1.20) , The Wannier center cross the reference 
line (red dashed line) once (odd times); (B) For the Normal 
insulating phase (A=-14.48, B=-18.0, C=-0.018, D=-0.594, 
M=2.766, A=1.20), The Wannier center cross the reference 
line (red dashed line) zero (even) times. 

To show these results more clearly, we glue 6 ~ — tt line 
and 6 = TT line together. Then the Wannier centers live 
on a cylinder surface. The corresponding results for the 
QSH phase are showed in Figl^JA). When moving from 
fcy = to TT we see that the two evolution lines of Wannier 
centers enclose the cylinder once and equivalently these 
evolution lines cross the reference line (the red dashed 
line) only once (odd times). By contrast, for the normal 
insulator phase, as shown in FiglJI^B), the two evolution 
lines of Wannier centers never cross the reference line. 
Therefore in the BHZ toy model for TI, the Z2 number 
calculated by our new method is consistent with the pre- 
vious conclusion. Next we will apply the method to more 
realistic models of insulating materials. 



B. CdTe and HgTe 

The CdTe and HgTe materials have a similar zinc- 
blende structure without bulk inversion symmetry. CdTe 



has an normal electronic structure, where the conduc- 
tion bands (Fg) have the s-like character and the valance 
bands have the p-like character(r8) through out the 
whole Brillouin Zone. In HgTe, the band structure is 
inverted in a small area near the F point, where the s- 
like Fg band sinks below the p-like Fg band. The band 
inversion at the F point changes the topological nature of 
the band structure and makes the HgTe to be the topo- 
logical insulator if a true energy gap is opened by the 
lattice distortion^(As pointed out in ref.? the uniaxial 
strain is applied along the [001] direction for HgTe by 
choosing the c/a ratio to be 0.98 and the energy gap is 
about 0.05eV at the F point). We then apply a tight- 
binding model^ to calculate the pattern of the Wannier 
center evolution 9 defined in Eq[5] and show the results in 
FigEl It is very clear that in the HgTe system, for k^ = 
the evolution line crosses the reference line (red dashed 
line) once (as shown in Figl^l^A)), while for fc^ = tt it 
never crosses (as shown in Figl31[B)). The above results 
indicate that in the case of HgTe the effective 2D systems 
for fixed kz = and tt are effectively 2D topological insu- 
lator and normal insulator respectively, which determines 
HgTe to be a strong 3D topological insulator—. A similar 
analysis can be also applied to CdTe and the results are 
shown in FigI3]JC) and (D). They clearly indicate that 
CdTe is a normal insulator. 



C. Bi2Se3 system 



Recently, the tetradymite semiconductors Bi2Te3, 
Bi2Se3, and Sb2Te3 have been theoretically predicted 
and experimentally observed to be topological insula- 
tors (TI) with a bulk band gap as large as 0.3eV in 
Bi2Se3^i^i23^i26^^ rpj^g gj_^gg^ surface state has been 

found by both ARFES^ii^ and STM2^, consistent with 
the theoretical results^S. 

Since the Bi2Se3 family has inversion symmetry, the Z2 
topological number can be easily calculated by the prod- 
uct of half the parities at each high symmetry points in 
the Brillouin Zonei^. Below we apply our new method to 
calculate the topological property of this system, using 
the tight binding model based on the Wannier functions 
obtained in reference^. We first perform the calculation 
for the Bi2Se3 without spin-orbit coupling: the results 
are shown in Fig|4lA,B). It is clear that the evolution 
lines never cross the reference line for both kz = and tt, 
indicating that the system is topologically trivial with- 
out spin-orbital coupling. When the realistic spin orbital 
coupling is turned on, as shown in Fig|4UC,D), the evo- 
lution lines cross the reference line once only in the case 
of A;^ = but not for kz ^ tt indicating the Bi2Se3 bulk 
material is a 3D strong topological insulator. 
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Figure 3: The evolution lines of Wannier centers for HgTe (A, 
B) and CdTe (C, D). For HgTe system, The Wannier center 
cross the reference line odd times in the kz = plane and 
even times in the = n plane, indicating HgTe is a strong 
topological insulator. For CdTe, the Wannier centers cross 
the reference line zero times for both kz — and n planes, 
indicating CdTe is a normal insulator. 




Figure 4: The evolution lines of Wannier centers for Bi2Se3 
system without(A, B) and with SOC (C, D). If we turn off 
SOC the system is a normal insulator and the Wannier centers 
never cross the reference line for both k- = and n planes 
as shown in (A) and (B), indicating it is a normal insulator. 
When the SOC is turned on the system is in a strong topo- 
logical phase, the Wannier centers cross the reference line odd 
times in /c, = and even times in k^ — it plane, 
as shown in (C) and (D), indicating the system is topologically 
non-trivial. 
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D. Bi2Te3 slab system 

As calculated by Liu et al^, upon reducing the thick- 
ness of Bi2Te3 and Bi2Se3 films, the topological nature 
of the system alternates between topologically trivial and 
non-trivial behavior as a function of the layer thickness. 
Liu et al. point out that the IQL Bi2Te3 slab is a triv- 
ial insulator and 2QL Bi2Te3 slab is a 2D topological 
insulator—. We apply our method to these systems. The 
evolution patterns for the IQL and 2QL Bi2Te3 slabs are 
obtained using the tight binding Hamiltonian developed 
in references^i^ and the results are summarized in Fig IS] 
In the IQL slab system (FiglHlIA)), the evolution pattern 
appears in a trivial manner while that of the 2QL slab 
system is non-trivial(Fig[S{B)). This is consistent with 
the conclusion based on the parity counting^. 



Sb shows clearly that it is in the normal insulator phase, 
which is also consistent with the parity counting. 





ky(27i) 



Figure 6: The evolution lines of Wannier centers for 2D single- 
bilayer Bi (A) and Sb (B) system, indicating the single-bilayer 
Bi is topologically nontrivial but Sb is topologically trivial. 
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Figure 5: The evolution lines of Wannier centers for IQL and 
2QL Bi2Te3 slab. (A) The IQL Bi2Te3 slab is in normal 
insulator phase. (B) The 2QL Bi2Te3 slab is in topological 
insulator phase. 



F. Graphene system 



E. Bi and Sb system 

Murakami pointed out that the Z2 topological num- 
ber is odd in the 2D bilayer bismuth system^. We ap- 
ply our method to this system using the tight-binding 
model developed in reference^, which faithfully repro- 
duces the bulk bismuth band structure. As shown in 
FiglSK, the band structure of bilayer bismuth is topolog- 
ically nontrivial, which is quite consistent with the previ- 
ous conclusion^. After that we apply the same method 
to calculate single-bilayer Sb, which has the similar lat- 
tice structure as bismuth, but with relatively weak SOC. 
As plotted in Figl5{B), the evolution pattern of bilayer 



Li 2005 Kane and Mele pointed out^ that there are 
two different phases in graphene, depending on the spin- 
orbital coupling \ji and the staggered sublattice potential 
parameter A^. The system is in a quantum spin Hall 
(QSH) phase when = O.li and normal insulator phase 
when A„ = 0.4i, where t is the hopping parameter. In 
the present study, we use the same parameters as Kane 
and Mele to calculate the corresponding Wannier center 
evolution pattern. The results are shown in FiglT] It 
can be easily found that there is one partner switching 
in Fig[71[A) but not in Fig[7]JB) indicating the former is 
topologically non trivial and the latter is trivial. 
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Figure 7: The evolution lines of Wannier center for graphene 
in the (A) QSH phase Au — O.lt and (B) the normal insulating 
phase Au = OAt. In both cases \so ~ 0.064 and Xr = 0.054. 

In conclusion, we have proposed a new equivalent 
expression for the Z2 topological invariance using the 
U{2N) non-Abelian Berry connection. Based on this new 
expression we calculated the evolution of the Wannier 
function center for several topological and normal insu- 
lating systems with or without inversion symmetry. We 
showed that for the nontrivial topological insulators, the 
Wannier function center have partner switching patterns, 
topologically different from the normal (trivial) insulat- 
ing systems. Additionally, we gave a proof that the new 
method is equivalent to the Z2 number proposed by Fu 
and Kane. 

Note During the preparation of this manuscript, we 
noticed the paper by Soluyanov and Vanderbilt^, where 
the construction of Wannier functions for Z2 topological 
insulators are discussed from a different view of point. 
We addressed in this paper that the real construction of 
Wannier functions is not necessary, while only the "Wan- 
nier representation" and corresponding Berry connection 
evaluated along the "Wilson loop" are essential keys in 
order to identify the topological nature. 
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Appendix A: Proof of gauge invariance of the D{ky) 
operator and its equivalence to the 17(2) wilson loop 

We first give the proof of gauge invariance of the D{ky) 
operator and its equivalence to the U{2) Wilson loop. 
The D{ky) operator is, at each ky^ the product: 

D{ky) = Fo,iFi^2F2^:iFN^-2,N^-iFN^-ifi (Al) 
where 

i^M^i = (™, ky\n, ky) (A2) 

kx^i = lAfc, Afc = and let kx^Nj, = to make k^ go 
back to the initial point. m,ri are running through from 
1 to the number of occupied bands, \n, k) is the n'th 
wavevector of energy En of the insulating Hamiltonian: 

H{k)\n,k) ^ En\n,k) (A3) 

For a fixed ky, the k^ = [— tt, tt] interval is broken up in 
Nx parts. This corresponds to periodic boundary condi- 
tions in the x direction and A^^; sites. D{ky) operator is 
then: 

Dmn{ky) = [Fo^iFi^2F2,3Fn^-2,N^-iFn^-1 ,o]mn 

= (m, fco|ni,fci)(ni,A:i|n2,fc2) ••• 

X {nN^-2, kN^-2\nN^-i,kN^-i){nN^-i,kN^-i\n, ko) 

(A4) 

where double index of rii implies summation. As the ky 
index is self-understood, we will suppress it at times and 
re-introduce it when necessary. 

Under a unitary transformation \n, k) become to : 

\n,k) -» M„t„(fc)|m, fc) 

and hence D[ky) matrix goes into 

F>mn{ky) M*^,^{ko){m' , ko\n[,ki) 

X M„.„,(fci) • • • AO,^_^„^_^(/cw^_2) 

X {n'N^^2^kN,-2\n'N^-.i,kN^~l)Mn'^^_^nr,^-iikN^-l) 
^ - 1 - 1 (^^^- - 1 ) ^^'k -l,kN^-lW, ko)Mn'n {ko ) 

=M^'miko){m',ko\ni,ki) ■ ■ ■ 

X {nN^-i,kN^-i\n',ko)Mn'niko) (A5) 

or 
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D{ky) ^ M\ko)D{ky)M{ko) 

which means that the eigenvalues of D{ky) (or its 
trace and determinant) remain unchanged under a 
U (occupied) gauge transformation. 
For infinitesimal Ak ^ 2tt, we have: 



-^i^+i ("^i kx.i, ky\n, kx,i+i, ky) 



A -1_ 4'"" At 



(A6) 



where = ('".fcx.^-fc^KI".^^^..''^)-!".^^...^^)) jg non- 

abelian U{occupied) gauge field. Hence we have that 
D{ky) is: 



^2 : ^1 



B{k2)Fl^,^B\k,) 



(B4) 



We now focus on the ky ^ or ky ~ n paths (say 
^ = 0), each of which has k^ going form — tt to tt, so 



that ky 



-ky, and compute the finite difference: 



7V,-1 



Diky = 0) = n F,,,+i 



i=0 

X F^AkfiFa^AkFAk,2Ak ■ 



■■F, 



3Ak,-2AkF-2Ak,-Ak 



-l)Afe, 



(B5) 



where again we are at fcy = for all the F's in the above, 
and the k in the above expresses the k^ coordinate. By 
the above F-k2.-ki = B{k2)F^^ ^.^B'^{ki), we have: 



D{ky) = 



Ai i^.\Ak 



(A7) 



This is just the {/(occupied) {not S'[/(occupied)) Wilson 
loop, where the contour Cfc^^is a contour at fixed ky which 
goes across the BZ in k^, i.e. goes from k^ = — tt to 
kx = TT, through kx = 0. 



Appendix B: Wilson loop and pfafHan topological 
invariant 

In this section, we will prove that the U{2) Wilson loop 
is related to the topological invariant of time-reversal 
topological insulators. In a system with time-reversal 
invariant we can relate the bands at k and —k through a 
unitary matrix B, 



\n,^k)^Bl^{k)T\m,k) 
with B{k) unitary and has the property: 

B{-k) = -B^{k) (B2) 
We have the following relation between F^^^k^ matrices: 



Afe,0 



B{Ak)Fl^kB\Q), 



F-2Ak,-Ak = S(2A/c)fI,. 2^,B^(A/c), 
Hence the Wilson loop above becomes: 



(B6) 



D = BC-^Ak)F^, 



Nx 



2-^k)^{IhL-i)Ak,^Ak 



St((_^_l)Afc), 



xB{2Ak)Fik^^^kB^Ak)BiAk)F^^^,B^iO)Fo,Ak 



X Fa 



fe,2Afe 

,Nx 



F2A 



fe,3Afe 



■F, 



(^-l)Afc,^Afc 



^(^^^)-F'(^_l)Afc,i:^Afe-^(V-2)Afc,(i:i^-l)Afc ' 



X FXk,2AkFo,Ak^'' i^)^0.AkFAk,2AkF2Ak,3Ak ' ' ' 



X F, 



-2}Ak,{i 



'l)Ak^(/ 



-l)Ak,^^Ak 



(B7) 



where we have used the fact that B'^{k)B{k) = I (unitary 
matrix) . We hence see that all the intermediate B matri- 
ces vanish with the exception of the ones at the inversion 
(Bl) symmetric points 0, ^f-Ak. This is obviously true for any 
time-reversal invariant contour. Moreover, it is sugges- 
tive that the two left-over matrices i?^(0), B{^Ak) be 
brought together, so we must commute 5^(0) all across 
the matrix chain. 

The matrix B{0) (and B{-k)) has the property that it 
is unitary and antisymmetric, per B{k) = B'^{—k) in 
ea(|B2|l . We then have: 



= {m',k2\fBm7n'B*^,f\n',ki) 



^mm'[-rki,k2) ^n'r 



(B3) 



BiO) 



e cr2 



(B8) 



The matrix -Ffc^fej, for k2 — fci <C tt, as show in eq (jA6[) 
has the following form: 



or equivalently: 



F^uk.=Srnn+AZ%{k2~h) 



(B9) 
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where k2 — ki = Ak. We decompose the U{2) gauge field 
into its abehan and non-abehan part; 



A 



ki.k2 — ^ki,k2 mn + ^kiM J"!" \^^^) 

1,2,3 and double index implies summation. 



where 

^fe/fe2' ^fei'^fe2^'* ^'^^ numbers. We have: 



y^kiMI - y^kiM) 



= <5„„„ + + Al-J>-\ar^^)[k, - k,) (Bll) 

We have that: 

iFlM)BHQ) = i?^(0)(/ + <« (fc2 - fci)/) 

+ All'i^^\ar{k2-k,)B\Q) (B12) 

we have commuted the first two terms easily as they are 
proportional to the identity, but the last term required a 
bit more work: 



.SU{2).i( i-.T 



T 



ay ay = -ayay 



a, ay 



(B13) 



and hence: 

AtT'(-rB\0) = e-^^4T-iar'^, 
Hence: 



(B14) 



{Fl,k2)BH^) 

= B\Q)[I + (<« / - 'V*)(fc2 - fci)] 

+ 2All%I{k2-k^)] (B15) 
We also have 



= / + {Ak,M)\k2 - fci) + '2A]!jll{k2 - k,)I 
« U + iAk.M)Hk2 - fci))(/ + 2Al^lik2 fci)) 



(B16) 



where in the limit of fc2 ^ fci ^ 27r (the case in all our 
terms) we neglect (fc2 — fci)^ order terms, and where 
Aki.k2 is the full C/(2) field strength. We have also used 
the fact that: 



AZ% ik2 - fci) = (m, fciKlTi, fc2) - \n, fci)) 

(m, fci|n, fc2) ~ 6mn 
= S„in - (JTI, fc2|n, fci) 
= -[(n, fci|m, fc2) - 6rnn]* 

= -(^rfeJ*(fc2-fci) (B17) 

For fc2 — fci ^ 27r, to give the above equation, we used 
the finite difference version of the equality: 

(to, k\n, k) ^ 5jnn {dk{m, k\)\n, k) + {m, k\{dk\n, fc)) = 
which reads by putting: 



and 



(9fci (to, fci|)|ri, fci) 



(to, ki\{dki |n, fci)) 



((m, fc2| - (to, fci|)|n, fci) 
fc2 — fci 

(to, fci|(|n, fc2) - |n, fci)) 
fc2 — fci 



(m, fc2|n, fci) + (to, fci|n, fc2) = 2(5,„„ 
After all this long detour, we have proved: 



We then return to the U{2) Wilson loop. 



(B18) 



(B19) 



D - Bi^Ak)F, 



pi 

'l)Ak,^Ak i 



-l)Ak 



-2)Ak,{- 

^ ■^(V-3)Afc,(^-2)Afc-^(^-4)Afe,(^-3)Afc 

^ -^Afe.2Afc F(^^AkBH0)Fo,AkFAk,2AkF2Ak ,3Ak ■ ■ ■ 

= B{^Ak)F^^_^^^^,^^i^F^^_^^^i^^^^_^^^^ 



^ -^('V-3)Afc,(^-2)Afc-^(^-4)Afe,(i^-3)Afc ' 



^Fl,,^2AkBH0)e^<'^'>'^''Fl^,Fo.AkFAk 

X i^(i^_2)Afc,(i^-l)Afc-^(i^-l)Afc,i^Afc 



2Ak ■ ■ ■ 



+ -+A 



)Ak 



X B{^Ak)B\0) 

where we have used C/n ai.'^i 
/, etc.... Hence: 



OAfc^O.Afc 



(B20) 

^Afc.2Afc^Afc,2Afc = 
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X B{tt)B^O) 



{Nx /2-l)Ak,{Nj: /2)Afc 



)Ak 



(B21) 



We note that the phase above is twice the C/(l) (ma- 
trix, not traced) phase picked up from to tt. Lets re- 
define it by re-expressing it from — tt to tt. i.e. the full 
abelian Berry for the interval considered. We have 



<W Afc = i ^(n, fciKIn, fc2> - \n, fci)) (B22) 

n 



-tr[F_k2.-ki - I] 
^tr[B{k,)Fl^,^B\h)-I] 



(B23) 



where we have used ea (jB4p . Since k2 — ki <^ 2tt, we can 
approximate B{k2) — B{ki) as small and write B{k2) = 
B{ki) + B{k2) - B{ki) to get: 



, Afc 

+ (B(fc2)-S(fci))F,^,fe,i?nfci)] 

- ^^'^[^A^.fe. - I + {B{k2) - B{k,))Fl^^^B\k,)] 

= \tr[Fk,M - I + {B{k2) - B{k,))Fl^^^B\k^)] 



A'^.l^lAk + -trmk2) - Bik,))Fl,.B\k,)] 

(B24) 



As B{k2) — B{ki) is considered small for k2 — fci <C 27r, 
we take 



{B{k2) - B{k,))Fl^^^ « i3(fc2) - S(fci) (B25) 

where we took Fj^^ w / if multiplied by another small 
number. Hence: 

Afc = A^lAk + ltrmk2) - B{k,))B\k,)] 

(B26) 

We then find: 



^V^O.Afc ^ ^Afc,2Afe ^2Afc,3Afe 

^-2A/c,-AA; + ^-Afc,0 + ^O.Afc ^ ^A/c,2Afc 



(N^/2-l)Ak,(N^/2)Ak 



dktr[B\k)V kB{k)] 



)Afc 



(B27) 



The first term is just the U{\) phase in the contour di- 
rection: 

A^'^^\k)dk. The Wilson loop is then: 

(B28) 

As B{k) is unitary, we also know that: 



tr[B'^{k)VkB{k)] = VklogdetB{k) 
and hence: 



(B29) 



The Wilson loop becomes: 



ldetB{0) 
detB{Tr) 
(B30) 



D = e 



-/in 



(B31) 



As we said before, B{0), B{tt) are unitary, 2 by 2 ma- 
trices, antisymmetric, so: 



B{0) = P/(i?(0)) 



B{7t) = Pf{B{Tl)) 



1 

-1 



1 

-1 



(B32) 



(B33) 



where Pf is the pfaffian of the matrix. We hence have: 



" idetBin) P/(B(0))^ ^^^^^ 

where I is the 2x2 identity matrix. 

We now make several observations. Obviously, the 
above equality is valid on both time-reversal invariant 
lines at fcy = and fcj, = tt, i.e. we can define two Wilson 
loops: 
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D{ky = 0) = e 



, detBjO.O) P/(B(7r,0)) r 
detB{n,0) P/(B(0,0))^ 



(B35) 



and 



Z?(fc, = tt) 



/ rfetB(0,7r) P/(i3(7r,7r)) ^ 
detB{TT,T!) P/(S(0,7r)) ^ 



(B36) 



Second, we notice that the U{1) phase factor is not just 
the usual abehan Berry phase but only half of it. Indeed, 
as per our definition: 



.m,n _ .(7(l),ir , 



SU{2)S, i 



(B37) 



This implies: 



\Y.{m,k\m,k) (B38) 



which has a 1/2 difference from the usual form. This 
difference is actually important. Define: 

<5{ky) ^ (j) dk^A^{k^, ky) ^logdct D{ky) (B39) 

We then have: 

/ ^kyHky)= VkJOgdct D{ky) 

Jo Jo 

= $(7r) - $(0) + 27riM„ (B40) 

where Af„ is the winding number of the phase <I>(7r). The 
phase is the sum of the phases 4>i{ky) and 4>2{ky) 

of the two eigenvalues of the Wilson loop both defined 
in the interval [0, 27r]. Each of these eigenvalues has a 
winding number which adds to M„, and the system will 
turn to be nontrivial if the system has a odd Af„. We 
now take the Wilson loop W from kx = — tt, it aX ky =^ 0, 
and then from k^ ~ tt, — tt at ky = tt: 

W = D{ky = Q){D{ky = tt))-! 



gi($(0)-$(A;)) X 



l detB{0,0) Pf{B{TT,0)) 
detBi7r,0) P/(B(0,0)) 



detB{Q,n) Pf{B{n,Ti)) ^ 
^ \] detB{Tr,Tr) Pf{B{0,Tr)) 

= pk(-l0^kylOgdctD{ky))+TTiM„ 



ldetB{0,0) Pf{B{TT,0)) 
detB{Tr,0) P/(B(0,0)) 



l detB{0,TT) P/(P(7r,7r)) 
deiP(7r,7r) P/(P(0,7r))' 



(B41) 



As such: 



P»(fcj, =0)e5(/frVfc„ iogdotD(fc,))p^^^ =7r))-i 



ldetB{0,0) P/(P(7r,0)) 
detB{TT,0) P/(P(0,0)) 



detB{0,TT) Pf{B{TT,7T)) 

y detB{TT,7r) Pf{B{0,TT)) 



(B42) 



We have proved that both -D(fci, = 0) and D{ky = tt) are 
proportional to unity matrix, up to a sign. For a smooth 
gauge, the difference in sign is taken by the contour term 

(B43) 

to give 



1 = 



/ detB{0, 0) P/(P(7r, 0)) / deiP(0, tt) P/(P(7r, tt)) 
Y detB{TT,0) P/(P(0,0)) y detB{n,TT) P/(P(0,7r)) 

(B44) 

which says that the pfaffian invariant is just the parity of 
the band switch number M„. For M„ odd, it is nontriv- 
ial. Note that although our proof above is explicit only 
for two occupied bands, it can be easily extended to the 
'^Noccupied band case when we realize that that case is 
just a tensor product (upon removing accidental degen- 
eracies) of Noccupied timc-revcrsal invariant multiplets for 
which the above expression applies. 



Appendix C: Wilson loop and the Z2 invariant 
expressed as an obstruction 

An alternative formulation of the Z2 invariant has been 
defined by Fu and Kanei, where the Z2 invariant is ex- 
pressed as an obstruction of the U{\) Berry's phase gauge 
field in half of the Brillouin zone. This approach has some 
convenience in its similarity with the Chern number for- 
mula of the quantum Hall states by Thouless et al^. The 
application of this approach to numerical calculation of 
the Z2 invariant in finite size systems has been studied by 
Fukui and Hatsugaii^ . Here we provide an alternative 
proof of the relation between our Wilson loop approach 
and the Z2 invariant through the obstruction formula. 
We start by reviewing the obstruction formulation of Fu 
and Kane^ (Appendix Al). Consider \n,k) as the occu- 
pied Bloch bands. We make the gauge choice 



\n,-k) ^T„rnT{\m,k)) 



(CI) 



with T an antisymmetric matrix satisfying = —1. 
Comparing to Eq. (|Bip . the gauge choice here corre- 
sponds to the requirement that Bnmik) is independent 
from k. More explicitly, with 2A^ occupied bands we 
can label the bands in pairs as \n,k) and \n,k) with 
n — 1,...,N. Time-reversal acts as \n,—k) = T\n,k), 



-k) 



-T\n,k). so that the wavefunctions in the 
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lower half Brillouin zone (BZ) defined by ky G [— tt, 0] is 
determined by those in the upper half BZ, denoted by 
Ti/2- With this gauge choice, for topological insulator it 
is not possible to make a continuous and single-valued 
choice of the wavefunctions in the whole Brillouin zone. 
However, it is always possible to define the wavefunctions 
continuously in the half BZ T1/2, so that all obstructions 
are pushed to the boundary between the two half BZs, 
i.e. the two lines ky = and ky = tt. In such a gauge 
choice, Fu and Kane shows that the Z2 invariant is given 
by an obstruction in the half Brillouin zone: 




A dl 



/ (fkF^y I mod 2 (C2) 



where Ai = — i ^jj(nfc|(9i|nfc) is the C/(l) part of the 
Berry phase connection. By first integrating over kx and 
define the U{1) Wilson loop 



/•TT 

^{ky)^(p dkxAx{kx,ky) 
The Z2 invariant can be expressed as 



(C3) 



(C4) 

If we don't have any restriction on the gauge choice 
(other than requiring that the wavefunctions are con- 
tinuous in Ti/2 so that $(fcj,) is continuous and well- 
defined), A can be 0, or any arbitrary integer. However, 
the gauge choice eg. ljCip removes this ambiguity. Con- 
sider the states \nk) for ky ~ 0. A gauge transformation 



\nk) 



nk) 



corresponds to a gauge transformation a 
which leads to the change in the flux 



(C5) 



^{ky=G)^^ky=Q) + 2N (h dkAvk^^o (C6) 



(a) 




0iC7r) + 27r 



Figure 8: (a) Definition of the iialf Brillouin zone and the 
paths to define Wilson loops, (b) Schematic picture of eigen- 
values (f)n{ky) for a two band model. 4ii{ky) has winding num- 
ber and (j)i{ky) has winding number 1, when <^i(0) and (j>i{n) 
are chosen to be in [0, 2tt) . 



However, to preserve the condition eq. (jCl[) we have to 
require 



(C7) 



so that 



dkxdxfk^fi = 2 / dxipk^,o = 2 {(p^fl - (po,o) (C8) 



Since we also have tp^^o s-nd ipQ^o = mod tt from ea. (|C7|l . 
the allowed change of ^{ky = 0) in the gauge transfor- 
mations that preserves the gauge choice ea. (|Cip can only 
be 



Hky 



0) ^ ^{ky 



0) + 4iV (v3^,o - VJo.o) 

0) + 4?™, n e Z (C9) 



The same is true for $(fcy = tt). Consequently, ^{ky = 0) 
and $(fcj, = tt) are well-defined modular 4tt, so that the 
Z2 quantity A defined by ea. (jC4|) is well-defined. 

Now we relate this result to the non-Abelian Wilson 
loop. When there are 2N bands occupied, a U{2N) Berry 
phase gauge field a"™ = —i{nk\di\mk) is defined. We can 
define the U{2N) Wilson loop along the same equal-fcy 
loops: 



W{ky) = Pe^fdk^<^^{k^^ky) e U{2N) 



(CIO) 



The U{1) gauge field is related to the U{2N) gauge field 
by Ai = Troi so that the U{\) flux ^{ky) is related to 
W{ky) as e**^'^"^ = detVF(fcy). Denote the eigenvalues 
of W{ky) as e^'t'-^'^^'y) with n = 1, .., 2A, we have $(fcy) = 
HiAniky) mod 27r. Thus 

^ = f dkydkJn{ky)-Y,{<i^n{T^)~Km] 

\ ri •'0 n / 

mod 2 (CU) 

Now we study the effect of the gauge choice eg. ljCip on 
the U{2N) gauge field. 



sT^ = -^(n,-fc|V 



= iTrnlTnpT ((pfc|Vk|/fc)) 



<) =(Ta^r-)_ (C12) 



W{-ky)^TW' {ky)T- 



(C13) 



Thus for ky = or tt, e**^"^'"'"-' is doubly degenerate. 
If we label a pair of degenerate eigenvalues by n and 
fi, the gauge choice ea. (|Cip corresponds to the choice of 
(f>n{ky) — (j)ff{ky). Indeed we see that if we make this 
choice, an ambiguity of 27r in leads to an ambiguity of 
Att in ^„ (/)„. Thus A defined in eq. (|Cll[) is well-defined. 
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To simplify the formula, we choose c/)„(0) and 4>n(j^) to 
be in [0,27r). Thus 




dkydkjniky) = 0„(^) ~ <^„(0) + 27rA/„ (C14) 



with Mn the winding number of phase 0„, which is 
equal to the number of times (/)„ crosses the line 0„ = 27r 
from below. For example in Fig[H] (b) (pi has winding 
number and 0j has winding number 1. In this way we 
get 



A = ^ Mn mod 2 (C15) 



The number M„ simply counts how many eigenval- 
ues 0„ crosses = 27r line (or any other reference line) 
from below. Thus the Z2 invariant is simply determined 
by the parity of the number of eigenvalue curves 4>n{ky) 
which crosses a reference line (j) = constant. 
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